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WILLMAN 1 - A PROBABLE DWARF GALAXY WITH AN IRREGULAR KINEMATIC DISTRIBUTION 
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ABSTRACT 

We investigate the kinematic properties and stellar population of the Galactic satellite Willman 
1 (Wil 1) by combining Keck/DEIMOS spectroscopy with KPNO mosaic camera imaging. Wil 1, 
also known as SDSS J1049+5103, is a nearby, ultra-low luminosity Milky Way companion. This 
object lies in a region of size-luminosity space (My ~ —2 mag, d ~ 38 kpc, n,alf ~ 20 pc) also 
occupied by the Galactic satellites Bootes II and Segue 1 and 2, but no other known old stellar 
system. We use kinematic and color-magnitude criteria to identify 45 stars as possible members of 
Wil 1. With a systemic velocity of ^hciio = —12.8 ± 1.0 km s _1 , Wil 1 stars have velocities similar to 
those of foreground Milky Way stars. Informed by Monte-Carlo simulations, we identify 5 of the 45 
candidate member stars as likely foreground contaminants, with a small number possibly remaining 
at faint apparent magnitudes. These contaminants could have mimicked a large velocity dispersion 
and abundance spread in previous work. We confirm a significant spread in the abundances of the 
likely Wil 1 red giant branch members ([Fe/H] = -1.73 ± 0.12 and -2.65 ± 0.12, [Ca/Fe] = -0.4 
± 0.18 and +0.13 ± 0.28). This spread supports the scenario that Wil 1 is an ultra-low luminosity 
dwarf galaxy rather than a star cluster. Wil l's innermost stars move with radial velocities offset by 
8 km s~ x from its outer stars and have a velocity dispersion consistent with km s~ 4 , suggesting 
that Wil 1 may not be in dynamical equilibrium. The combination of the foreground contamination 
and unusual kinematic distribution make it difficult to robustly determine the dark matter mass of 
Wil 1. As a result, X-ray or gamma-ray observations of Wil 1 that attempt to constrain models of 
particle dark matter using an equilibrium mass model are strongly affected by the systematics in the 
observations presented here. We conclude that, despite the unusual features in the Wil 1 kinematic 
distribution, evidence indicates that this object is, or at least once was, a dwarf galaxy. 

Subject headings: galaxies: star clusters — galaxies: dwarf — galaxies: kinematics and dynamics — 
galaxies: individual (Willman 1) 



1. INTRODUCTION 

Since 2004, over a dozen Milky Way satellites 
have been discovered via slight statistical overdensi- 
ties of individual stars in the Sloan Digital Sky Sur- 
vey (SDSS) catalog and confirmed b y both fo llow- 
up im a ging and spectrosco py (e.g. Willman ~et all 
2QQ5 JE1 IZucker et all l2006aTa IBelokurov et alT 2006, 
2007t iSakamoto k Hasegawa]l2006|: llrwin et all 12007 : 
Walsh et all 120071: IBelokurov et all 12008. 200^ These 
satellites are dominated by old stellar populations and 
have absolute magnitudes of —8 < My < —1 mag. 
Their median My is ~ —4, less luminous than the me- 
dian observ ed for Milky Way globular clusters (GCs; 
IHarrisl[l996[ ). Stellar kinematics consistent with mass- 
to-light (M/L) ratios > 100 demonstrate that most of 
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these objects are dark matter dominated dwarf galaxies 
(Muhoz et al. 20061 iSimon fc Gehal [20071: iMartin et all 
2007; Strigari et al. 2008a). 

Four of the new Milky Way companions - Willman 1, 
Bootes II, Segue 1 and Segue 2 - contain L < lOOOi© and 
have been difficult to classify. With estimated rhaif of 20 
- 40 pc, these four objects lie in a gap between the sizes 
of known old stellar populations (Milky Way GCs and 
dwarf spheroidals) in size-luminosity space. They are 
less luminous than all but three known objects classified 
as globular cluste rs, providing few stars bright enough for 
kinematic study (iWillman et al.ll2005al IBelokurov et all 
[20071 iWalsh et alX 2008; B elokurov et all 120091 ) . More- 
over, their proximity to the Milky Way (d <; 40 

kpc) and their possible embedding in the Sagittarius 
strea m (Bootes II and Segue 1 and 2, IBelokurov et all 
20091: iNiederste-Ostholt et all 120091 - although see 



Law fc Ma icwski 2010) complicate the interpretation of 



their observed properties. 

Measuring the dark mass content of satellites with 
My > —3 is a critical ingredient to our understand- 
ing of the size and mass scale of dark matter cluster- 
ing, the abundance and distribution of dark matter ha- 
los, and the ex t reme l ow m ass limit of ga laxy formation. 
lKoposovet"aTI (|2007l ) and IWalsh et all (120091) showed 
that Milky Way companions fainter than My ~ — 3 could 
not have been discovered at all in SDSS if they are more 
distant than ~ 50 kpc from the Sun. They may thus rep- 
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resent the tip of an iceberg of such o bjects around the 
Milky Way (e.g. iTollerud et all 120081) . Moreover, these 
objects have been shown to be excellent targets for obser- 
vations seeki ng the gamma-ray sig nature of annihilating 
dark matter (jStrigari et al.ll2008btlGeha et al.ll2009f ). 

Two primary lines of reasoning have been used to argue 
for dark matter content in, and thus a galaxy classifica- 
tion for, these four extreme objects: i) a kinematic dis- 
tributio n that is unbound without dark matter (Segue 1 
and 2: iGeha et al.ll200a iBelokurov et al.ll2008| ). and/or 
ii) a spectroscopically observed spread in [Fe/H] , which is 
not expected in purely stellar systems with a total mass 
as low as the observable baryonic masses o f the ultra- 
faint dwarfs (Willman 1: IMartin et aLll2007D . Thus far, 
the strongest evidence for a substantial dark mass com- 
ponent has been provided by t he line-of-sight velo cities 
of stars in the Segue 1 object. IGeha et al.l (120091 ) ana- 
lyzed 24 member stars observed with DEIMOS to find 
(M/iy) contra i > 1000. Simon et al. (in preparation) 
confirm this result with a larger sample of 71 member 
stars. Circumstantial lines of reasoning have also been 
used to argue that Boot es II may contain a substantial 
dark matter component (jWalsh et al.ll2008J ). 

The reliability of these kinematic or spectroscopic 
[Fc/H] analyses of nearby dwarf galaxies hinge on having 
samples of member stars that are as contamination free 
as possible, and a quantitative calculation of the unavoid- 
able contamination that may be present. Contaminants 
may be stars from the Milky Way, from unrelated stellar 
streams such as Sagittarius, or from an unbound com- 
ponent of the dwarf galaxy itself. With a set of only 
^10-50 member star velocities, a small number of 
interlopers could artificially inflate the observed veloc- 
ity dispersion, leading to an overestimate of the mass- 
to-light ratio. With only a few stars in each of these 
systems bright enough for a spectroscopic [Fe/H] mea- 
surement, just one or two interloper stars at bright ap- 
parent magnitudes could mimic an [Fe/H] spread. Fore- 
ground Milky Way thick disk and halo stars (at the pho- 
tometric depths reachable by spectroscopy) contaminate 
the color-magnitude diagrams (CMDs) of these extreme 
satellites. Segue 1 and Bo otes II have syste mic veloci- 
ties of -206 ± 1.3 km s^lfGeha et al.ll2009h and -117 
km s~ 1 (|Koch et al.ll2009f ). respectively, which are offset 
from the majority of thick disk and halo stars. How- 
ever, both Willman 1 (v. 



.Sl/S 



= -13.3 ± 2.5 km s" 



Martin et all 120071) a nd Segue 2 (v sys = -40 km s -1 , 
Bclokurov et alll2009D have systemic velocities that sub- 



stantially overlap with the velocities of Milky Way stars, 
making the identification of interlopers particularly dif- 
ficult. 

Willman 1 (Wil 1; SDSS J1049+5103), located at 
(o^ooo, ^2000) = (162.343°, 51.051°), is an old, metal- 
poor Milky Way sa tellite at a distance of 38 ± 7kpc 
with M v ~ -2 mag ([Willman et al.H2005aHMartin et all 
2008a). Based on equilibrium models of its mass, this 
object h as been claimed t o have a high d ark matter 
content (IGeha et all 12001 IWolf et all l2010l) . A high 
dark matter content plus its relative proximity would 
make it a promising s ource of gamma-rays from anni- 
hilating dark matter (jStrigari et al.l l2~008b[ ). As a re- 
sult of this prediction, several groups have attempted 
to investigate the particle nature of dark matter via 



gam ma-rav (lEssig et all 120091: lAliu et all l2009D and X- 
ray ([Loewenstein fc Kusenkol 120101 ) observations. How- 
ever, its possibly irregular spatial distribution supports 
the idea that it is t idally disturbed (|Willman et alll2006l : 
IMartin et~a l. 2007) which could mean that its kinematics 
do not faithfully tra c e its g ravitational potential. Al- 
though IMartin et all (|2007f ) argued that Wil l's clas- 
sification of a dwarf galaxy was supported by an ob- 
served spread in the [Fe/H] of its member stars, such 
spectroscopic studies of Wil 1 suffer particularly from 
the presence of interloper stars. Its systemic velocity is 
similar to the velocities of Milky Way foreground stars 
with col ors and magnitudes consistent with Wil 1 mem- 
bership. Sieg el etall (|2008l) showed that c ontamination 
from M ilky Way stars was a problem in the IMartin et all 
(2007) study and found that 2 - 5 of the 7 Wil 1 spectro- 
scopic red giant branch "members" were actually Milky 
Way foreground stars. These foreground stars might 
have generated the apparent spread in [Fe/H]. 

To address the present uncertainties in the nature of 
the Wil 1 object, we present DEIMOS observations of 
45 probable member stars. We carefully characterize the 
possible contamination in this sample, and then use it to 
study the abundances and kinematics of stars in Wil 1. 
In §[2] we discuss target selection and data reduction for 
our DEIMOS slitmasks. Selecting Wil 1 member stars, 
including a detailed discussion of foreground contamina- 
tion is presented in §[3j In §[4j we discuss whether Wil 1 
appears to be a star cluster or a dwarf galaxy, and then 
analyze Wil l's kinematics in §[5] 

2. DATA 
2.1. Photometry 

The data are from wide-field imaging of Wil 1 with the 
mosaic imager on the 4m at Kitt Peak National Obser- 
vatory (KPNO) on April 7 and 8, 2005. The reduction, 
photometry and calibration of these da ta used in this 
paper is identical to that presented in Wil lman et all 
(2006). We briefly describe here the reduction, pho- 
tometry, and p hotometric calibration of these data. See 
Will man et al.l (|2006[ ) for more details. 

Ten 600s exposures were taken through each of the 
SDSS g and r filters, with seeing varying between 1.2" 
and 1.4". Photometry was performed on each individual 
exposure and then averaged to yield a stellar catalog with 
excellent photometry and star-galaxy separation at the 
apparent magnitudes we are investigating. 

Stellar magnitudes were photo metrically calibrated 
with the SDSS Data Release 4 (|Adelman-McCarthvl 
I2006D stellar catalog. The apparent ma gnitudes were 
then corrected for extinction using the Schlegc l et all 
(1998) dust map values given in the SDSS catalog. 
The average E(g-r) along the line of sight to Wil 1 is 
0.014. All magnitudes in this paper are de-reddened; 
we use the subscript "0" to denote that the colors and 
magnitudes are de-reddened. At magnitudes brighter 
than ro = 17.0, we replaced the KPNO photometry 
with the PSF photometry from SPSS Data Release 7 
(|Abazaiian fc Sloan Digital Sky Survevll2009l ). 

Figure [1] shows the CMD of all point sources within two 
elliptical half-light radii of the ce nter of Wil 1 , using th e 
structural parameters derived by IMartin et all (|2008bl) . 
The points in this CMD delineate the main sequence of 
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Fig. 1. — Dereddened color-magnitude diagram of all stars within 
two elliptical half-light radii of the center of Willman 1 from KPNO 
g- and r-band photometry . We used (position angle, ellipticity, 
f'half) = (77,0.47,2.3') from Mart in et all ( 2008a) to calculate half- 
light distances. The region inside the dotted boxes is the location 
of our highest priority spectroscopic selection criteria, hereafter 
referred to as the color criteria used to identify stars possibly be- 
longing to Wil 1. The sizes of color and magnitude uncertainties 
are shown by the crosses on the left of the CMD. 

Wil 1, with candidate red giant branch (RGB) stars vis- 
ible. Dotted boxes outline the liberal color-magnitude 
selection that we will use in the rest of this paper as the 
color-magnitude requirements for possible Wil 1 mem- 
bership. 

We applied the color-magnitude selection shown in 
Figure 1 to this photometric catalog to calculate a re- 
vised center of Willman 1. We began with the cen- 
ter, bas ed on the much shallo wer SDSS dataset, calcu- 
lated by iMartin et al.l (|2008bl ). and then iteratively cal- 
culated the average position of stars within 2 arcmin- 
utes of the center until we converged on (a2000: ^2000) = 
(162.3397,51.0508). We will use this center for the rest 
of the paper. 

2.2. Spectroscopic Target Selection 

Stars in Wil 1 were targeted for spectroscopy using 
the photometric catalog described in the previous sec- 
tion. We set the target priorities to preferentially ob- 
serve stars with a high likelihood of being Wil 1 mem- 
bers based on their color, magnitude and spatial position. 
First priority was given to stars that (1) spatially over- 
lap the main body of Wil 1 and (2) reside within regions 
of the color-magnitude diagram that are consistent with 
the Main Sequence (MS) and turnoff, horizontal branch, 
and red giant branch of an old stellar population at the 
distance of Wil 1. These color-magnitude criteria are 
shown by the dotted lines overplotted on Figure [TJ We 
chose to implement liberal, rectangular color-magnitude 
criteria to include Wil 1 member stars with a range of 
possible [Fe/H] and ages in our spectroscopic sample. 
Second priority was given to stars occupying a similar 
color-magnitude region, independent of spatial location. 
All remaining stars were assigned third priority. Within 
each of these three tiers, stars were further prioritized by 
their apparent magnitude, with the brightest stars re- 



ceiving highest priority. An average of 100 slitlets were 
placed on each mask (see Table 1). 

2.3. Spectroscopy and Data Reduction 

Four multislit masks were observed for Willman 1 us- 
ing the K eck II 10-m telesc ope and the DEIMOS spec- 
trograph (|Faber et al.l l2003). Three masks were observed 
on the nights of November 20-22, 2006, the fourth was 
observed on March 20, 2007. Exposure times, mask po- 
sitions and additional observing details are given in Ta- 
ble 1. The masks were observed with the 1200 line mm -1 
grating covering a wavelength region 6400 — 9100A. The 
spatial scale is 0.12" per pixel, the spectral dispersion of 
this setup is 0.33A, and the resulting spectral resolution 
is 1.37A (FWHM). Slitlets were 0.7" wide. The seeing 
conditions during both runs were on average ~ 0.75". 
Despite the similar observing conditions, few spectra 
were usable from the fourth mask because the targeted 
stars were fainter. The minimum slit length was 4" to al- 
low adequate sky subtraction; the minimum spatial sep- 
aration between slit ends was 0.4" (three pixels). 

Spectra were reduced using a modified version of the 
spec2d software pipeline (version 1.1.4) developed by the 
DEEP2 team at the University of California-Berkeley 
for that survey. A detailed descripti on of the two- 
dimen sional reductions can be found in iSimon fc Gehal 
(2007). The final one-dimensional spectra are re- 
binned into logarithmically spaced wavelength bins with 
15 km s" 1 per pixel. 

2.4. Radial Velocities and Error Estimates 

We measure radial velocities a nd estimate velocity er - 
rors using the method detailed in ISimon fc Gehal (|2007f ). 
We refer the reader to this paper for a description of the 
method and only highlight the important steps below. 

Radial velocities were measured by cross-correlating 
the observed science spectra with a series of high signal- 
to-noise stellar templates. The templates were observed 
with Keck/DEIMOS using the same setup as described 
in § 12.31 and cover a wide range of stellar types (F8 to 
M8 giants, subgiants and dwarf stars) and metallicities 
([Fe/H] = -2.12 to +0.11). We calculate and apply a 
telluric correction to each science spectrum by cross cor- 
relating a hot stellar template with the night sky ab - 
sorption lines following the method in lSohn et al.l (J2007). 
The telluric correction accounts for the velocity error due 
to mis-centering the star within the 0.7" slit caused by 
small mask rotations or astrometric errors. We apply 
both a telluric and heliocentric correction to all veloci- 
ties presented in this paper. 

It is crucial to accurately assess our velocity errors be- 
cause the internal velocity dispersion of Willman 1 is 
expected to be comparable to the DEIMOS velocity er- 
rors associated with individual measurements. We de- 
termine the random component of our velocity errors us- 
ing a Monte-Carlo bootstrap method. Noise is added 
to each pixel in the one-dimensional science spectrum. 
We then recalculate the velocity and telluric correction 
for 1000 noise realizations. Error bars are defined as 
the square root of the variance in the recovered mean 
velocity in the Monte- Carlo simulations. The system- 
atic contribution to th e velocity error was determined by 
ISimon fc~G cha ( 2007|) to be 2.2 km s" 1 based on repeated 
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independent measurements of individual stars, and has 
been confirmed by a larger sample of repeated measure- 
ments. While we did not place stars on multiple masks in 
the Willman 1 dataset, the systematic error contribution 
is expected to be constant as this is the same spectro- 
graph setup and velocity cross-correlation routines are 
identical. We add the random and systematic errors in 
quadrature to arrive at the final velocity error for each 
science measurement. 

The fitted velocities were visually inspected to ensure 
the reliability of the fit. Radial velocities for 111 stars 
of the 423 objects targeted passed this visual inspection. 
For the rest of this paper, we only consider the 97 of those 
111 stars that have spectra with S/N > 2, and that have 
velocity uncertainties of less than 10 km s . The me- 
dian velocity uncertainty of these 97 stars is 3.5 km s _1 . 
The median velocity uncertainty in the subsample of 
those stars most likely to be Wil 1 members is 4.7km s _1 , 
because the faintest of these 97 stars are dominated by 
main sequence stars belonging to Wil 1. All velocity 
histograms shown in this paper thus have a binsize of 
4.7km s -1 . 

2.5. Comparing relative and deriving absolute 
abundances of individual stars 

In §3.3 and §3.4 we will use the relative stellar abun- 
dances imprinted on these spectra to identify likely bright 
member stars, and to compare with stars observed by the 
SEGUE survey. In §4.1 we will then calculate the abso- 
lute [Fe/H] abundances of three bright member stars. 
We describe here the techniques used to measure those 
relative and absolute abundances. 

To assess relative stellar abundances, we use the re- 
duced equivalent width (W) of the Ca II lines. We 
utilize two functional forms of W, one whic h has been 
calibr ated for low meta llicity stars (Starkc nburg et al.l 
I2010D and the standard IRutledge et all (|1997t ) calibra- 
tion. For both, we measure the CaT lines at 8498, 8542, 
and 8662 A using the cont i nuum and line definitions de- 
scribed in IRutledge et al.l J1997) to calculate ECa. We 
determine the uncertainty on ECa with the Monte Carlo 
method described above. Added in quadrature to the 
Monte Carlo uncertainties is a systematic uncertainty of 
0.25 A, determined us ing the same method described in 
I Simon fc Gehal (|2007l ). but using a larger sample of re- 
peated measurements. 

Throug hout the paper, we primarily report the value 
of the IRutledge et al.l (|1997 ) definition of W' such that 
ECa = 0.5EW +EW +0.GEW and W = 

8498A 8542A 8662A 

ECa— 0.64(Vhb — V). V is the V-band magnitude of each 
RGB star and Vhb is the magnitude of the horizontal 
branch. To obtain V , we convert our SDSS g- and r-band 
magnitudes into V-band using the ph otometric transfor- 
matio ns given in Tables 1 and 2 of iBlanton fc Roweli 
(|2007[) . The apparent V magnitude of the one spec- 
troscopically confirmed star in the flat part of Wil l's 
horizontal branch is 18.45 mag. Although this approach 
minimizes assumptions, using a single star to determine 
Vhb introduces a possible offset into our W calculations; 
if the star is variable then its current apparen t magnitude 
may n ot equal its average value. However, iSiegel et al.l 
(2008) found no RR Lyrae stars in Wil 1 in their time- 
series imaging of the object. Regardless, a small shift the 



Vhb we use to calculate W would not affect any conclu- 
sion of this paper. 

We also use th e met allicity calibration detailed in 
IStarkenburg et al.l (|2010D to calculate [Fe/H] from W, 
where W is defined as 0.43(V - V aB ) + SCa (2+ 3) - 
2(ECa^ 2 1 + 5 3) ) + 0.034(U - Vhb), and ECa (2+3) = 
EW • +EW °. This study excluded the Ca 

8542A 8662A J 

II line at 8498A because it is the weakest of the 
three lines and has been shown to contribute more 
to the relative uncertainty tha n the two stronger lines 
(| Armandroff fc Da Costal [19911) . The CaT- [Fe/H] re- 
lation presented in IStar kenburg et al.l (HpiO) was cali- 
brated for RGB stars in the metallicity range —4.0 < 

[Fe/H] 0.5 and which lie within -3 < {V - Vhb) < 

or -3 < (M v ) < 0.8. 

To derive the absolute abundances of bright mem- 
ber stars, w e use the spectral synthesis method of 
IKirbv et all (|2008aft (KGS08). This method relies 
on comparing an observed medium-resolution spectrum 
with a grid of synthetic spectra covering a range of effec- 
tive temperature (T e ff), surface gravity (log <?), and com- 
position. All of our spectra with high enough S/N for this 
technique are bright enough to be in the SDSS database, 
so we determine V and I for each star by transforming 
th e SDSS qri magni tudes with the relationship derived 
bv lJordi et all (|2006[ ). KGS08 found that uncertainty in 
the measured colors, ages, and alpha- abundances of stars 
does not substantially affect the estimated T g and log 
g. The best-matching composition is then found by min- 
imizing residuals between the observed spectrum and a 
smoothed synthetic spectrum matched to the DEIMOS 
spectral resolution. To derive error bars on the best- 
matching [Fe/H] , we calculate % 2 contours for every star 
by allowing both [Fc/H] and T e ff to vary. We separately 
varied T off by ±125 K and 250 K and log g by ±0.3 dex 
and 0.6 dex. The surface gravity makes almost no differ- 
ence in [Fe/H] because there are no visible ionized Fe lines 
in red giants in the observed spectral range. Changing 
T eff by ±125 K yields a <J[Fe/H] of 0.13 dex and by an un- 
realistically large ±250 K yields a J[Fe/H] of 0.26 dex for 
each star we study in this paper. KGS08 found that their 
technique measured [Fe/H] with 0.25 (0.5) dex accuracy 
on spectra of Galactic globular cluster stars with S/N ~ 
20 (~ 10) A -1 . The success of this method has been con- 
firmed by a comparison with high-resolution abundances 
of over 100 stars in GCs, the halo field, dwarf galaxies, 
and of six RGB st ars with -3 3 < [Fe/H] < -2.3 in the 
ultra-faint dSphs (IKirbv et al.ll2T)0Ml2l)irl . 

3. A SPECTROSCOPIC SAMPLE OF WIL 1 STARS 

To study the stellar population and kinematic prop- 
erties of Wil 1, we need to identify a sample of mem- 
ber stars with minimal contamination from interlop- 
ers. Wil 1 lies at relatively high Galactic latitude (l,b) 
= (158.6°, 56.8°). However, its system ic velocity is ~ 
-12.3 ±2.5 km s _1 ([Martin et al.ll2007D and the median 
velocity of all Milky Way stars in the direction of Wil 
1 is —15.0 km s" 1 using the Besancon Galaxy model 
(| Robin et al. 2003). In addition to overlapping in veloc- 
ity, stars in Wil 1 have colors and magnitudes very similar 
to the colors and magnitudes of stars in the Milky Way's 
thick disk and halo. 

In this section, we use color, magnitude, and velocity 
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to select an initial sample of 45 Wil 1 candidate member 
stars. We then derive the predicted contamination from 
Milky Way foreground stars within this member sample. 
Finally, we identify 5 stars that are probable interlopers 
and present the properties of the remaining 40 probable 
Wil 1 member stars. 

3.1. A Color-Magnitude-Velocity Sample of Wil 1 
Candidate Stars 

The first required criterion for Wil 1 membership is 
having a color and magnitude consistent with the stel- 
lar population of Wil 1. We use a loose color- magnitude 
(CM) selection, as shown by boxes overplotted with dot- 
ted lines on Figure [TJ We used this loose cut, rather than 
proximity to a fiducial cluster sequence, to avoid making 
a priori assumptions about the stellar population of Wil 
1. 58 of the 97 stars in our spectroscopic catalog satisfy 
these CM criteria. 

The second required criterion for Wil 1 membership 
is a velocity consistent with belonging to Wil 1. Fig- 
ure [5] shows the heliocentric velocity histogram of these 
58 stars, along with the velocity histogram of the 39 
stars that do not satisfy the color-magnitude selection 
for Wil 1 membership. The velocity distribution of the 
58 CM selected stars is strongly peaked, with 45 stars 
between —30 and km s _1 . We identify these 45 color- 
magnitude-velocity (CM-V) selected stars as likely Wil 
1 members. This does not necessarily mean that none 
of the 13 CM selected stars with outlying velocities are 
physically associated with Wil 1. However, the spatial 
distribution of those 13 stars at outlying velocities is not 
clustered around the Wil 1 center. 

We present in Table 2 the equatorial coordinates, r 
magnitudes, g—r colors, heliocentric velocities, and spec- 
tral S/N of the 45 CM-V selected Wil 1 member stars. 
We also include the CaT W (and uncertainty) for the 
15 possible red giant branch, as calculated in § 12.51 Ta- 
ble 3 contains the same data (but not W) for the 52 
non-member stars. 

3.2. Predicting the Number of Interlopers in the 
Color- Magnitude- Velocity Sample 

Figure [3] shows a CMD of the stars in our spectroscopic 
catalog. Filled symbols represent the 45 candidate Wil 
1 members selected in § 13.11 and open symbols represent 
the 52 foreground Milky Way stars. The number of open 
symbols overlapping with the filled symbols shows that 
shows that 40% of stars with colors and magnitudes con- 
sistent with the red giant branch of Wil 1 are foreground 
stars belonging to the Milky Way. These foreground stars 
were only identified because their line-of-sight velocities 
were different than those of Wil 1 stars. The median 
velocity of Milky Way stars passing the CM criterion 
for membership is —35.7 km s _1 (based on the Besancon 
Galaxy model), with 16% of these having —30 < v\ os < 
km s _1 . How many Milky Way interlopers remain in the 
CM-V sample of 45 candidate Wil 1 members? 

We simulate the number of interloper stars expected 
among the 45 candidate members using the Besancon 
Galaxy model. Because photometric stu dies suggest the 
presence of tidal feature s around Wil 1 (Wil lman et al.l 
120061: iMartinet aTll2007ft . we first predict the number of 
Milky Way contaminant stars without assuming that all 
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Fig. 2. — Velocity distributions of: the 58 stars that satisfy our 
Wil 1 color-magnitude selection criteria (open) and the 39 stars 
that do not satisfy these criteria (grey filled) . The dotted lines show 
the velocity range of — 30 < fhelio < km s used to select Wil 1 
member stars. Binsize is 4.7 km s" 1 , the median velocity error of 
the 58 stars passing the color-magnitude criteria for membership. 
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Fig. 3. — Color-magnitude diagram of the 97 stars with 
DEIMOS/Keck velocities. Open symbols show Milky Way stars. 
Filled symbols show probable Wil 1 member stars, as selected by 
color-magnitude and velocity (—30 < V < km s -1 ) criteria. Tri- 
angles, circles, and squares highlight stars belonging to the bRGB, 
fRGB, and MS/BHB sub-samples used to characterize foreground 
contamination. 5-point stars show those stars that did not satisfy 
the initial color-magnitude cut for membership. 

CM selected stars outside the Wil 1 velocity peak belong 
to the Milky Way. We instead use the Besancon model to 
predict the absolute number density of Milky Way stars 
satisfying the color-magnitude-velocity criteria for can- 
didate members. The predicted number of contaminant 
stars thus rests on the assumptions that the velocity dis- 
tribution of Besancon model stars and the absolute num- 
bers of stars in the Besancon model are correct. We later 
verify that this yields a reasonable prediction. 
The primary ingredients in our calculation are: 
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1- "-fg,vci, the projected number density of Milky Way 
stars in the Besancon model satisfying the CM- 
V criteria for Wil 1 membership. We calculated 
^fg,vei and its dispersion in 1000 small helds ran- 
domly placed in a 1 square degree Besancon simu- 
lation centered on the position of Wil 1. To do this, 
we shuffled the RAs and Decs of Besancon model 
stars before selecting each random field. The ran- 
dom helds each had an area approximately equal 
to that of our spectroscopic survey footprint. Be- 
cause the CM cuts applied to our data were liberal, 
we simply used the model CFHT-Megacam g and 
r magnitudes as a proxy for the observed SDSS 
g and r magnitudes. We convolved 4.7 km s _1 
measurement uncertainties, the median for the 45 
candidate members, to the model velocities of each 
Besancon star. The average number of possible in- 
terlopers in the CM-V sample within a given area 
of sky, A, is then iV cont:VO i = A * nf gtVe i- 

2. /targ, the fraction of stars in our photometric cat- 
alog satisfying the CM criteria for Wil 1 member- 
ship that also end up in our spectroscopic catalog 
of 97 stars. Not all stars satisfying the CM criteria 
for membership were targeted, and not all targeted 
stars had spectra with high enough S/N to be in 
the final spectroscopic catalog. To derive the true 
number of Milky Way interlopers, iV contiVO i needs to 
be multiplied by /targ- Because both the density of 
stars and fraction of sky covered by observations 
decreases with increasing Wil 1 distance, /targ is a 
function of distance from the center of Wil 1. We 
thus calculate /targ in each of three distance ranges: 
0-2 r ha if, 2-4 r ha if and 4-6 r ha it . 

The contamination in the Wil 1 candidate member 
sample is expected to be primarily composed of stars be- 
longing to the Milky Way's thick disk and halo. Because 
the relative number of thick disk and halo stars varies 
across the CMD (thick disk stars dominate at brighter 
apparent magnitudes, halo stars at fainter) and because 
these two galaxy components have different velocity dis- 
tributions, nfg.voi is a function of apparent magnitude. 
Because stars were prioritized by apparent magnitude in 
target selection, /targ is also a function of apparent mag- 
nitude. We therefore separately estimate the foreground 
contamination in three subsets of the Wil 1 population: 
i) the five relatively bright red giant branch (bRGB; ro 
< 19.5 mag), ii) the 10 faint red giant branch (fRGB; 
19.5 < 7'o < 21.3 mag); and hi) the 30 main sequence 
(MS) and blue horizontal branch (BHB) stars. 

The left panel of Figure 2] shows /targ as a function of 
distance from Wil 1 for the bRGB, fRGB, and MS/BHB 
subsamples. The overall target efficiency at all distances 
is much lower than one, because many stars at the faint 
end of our acceptable magnitude range (r = 23) for tar- 
get selection were not measured with high enough S/N 
for a robust velocity to be extracted. Many stars remain 
to be observed in the center of Wil 1 if such faint mag- 
nitudes can be reached. The target efficiency for fRGB 
stars remains high out to large distances from Wil 1, even 
though our fractional spatial coverage beyond 3rh a it is 
quite low. This is because our mask placement included 
two out of the six total stars between 4 and 6ri la if that 



are consistent with the fRGB of Wil 1. Although target 
efficiency is low at large d, only 4% of stars are expected 
to lie beyond 3 elliptical half-light radii of a system well- 
described by an exponential profile (see also § 14.21) . 

The average expected number of contaminant stars in 
each of the bRGB, fRGB, and MS/BHB CM-V subsam- 
ples in each annulus is: 

-^*cont,vcl — /targ X ^-annulus X ^fg T vcli \*-) 

and the average dispersion in this number is: 



^cont.vcl — Jtarg X ^iannulus X (Xfg : vcb 



(2) 



where A an nuius is the area of each elliptical annulus, 
and /targ i rifg V ei and <7f g vc i are calculated separately for 
each of the bRGB, fRGB, and MS/BHB subpopulations. 
We predict a fraction of one interloper star within each of 
the three subpopulations within each of the three annuli. 
The probability that any individual star is an interloper 
increases with distance and is shown in the bottom left 
panel of Figure |3J 

We sum the predicted average numbers (and disper- 
sions, in quadrature) of interlopers in each annulus, to 
derive the average total numbers (and dispersions) of in- 
terloper stars expected in each of the bRGB, fRGB, and 
MS/BHB subsamples and find 1.5 ± 1.0, 0.7 ± 0.5, and 
0.6 ± 0.3 stars, respectively. Although the number of 
stars in our sample increases at fainter apparent mag- 
nitudes, the number of expected contaminants is low in 
the 30 possible main-sequence members. This small pre- 
dicted contamination results from the low overall target 
efficiency combined with the large number of Wil 1 stars 
in its main sequence relative to its bRGB and fRGB. To 
convert these fractional numbers of stars into a physical 
number of stars that may be in our sample, we used the 
IDL function POIDEV to generate 10 5 Poisson random 
deviates with the predicted mean and dispersion in the 
number of interloper stars. The result of this simula- 
tion is shown in the right panel of Figure SJ 50% of all 
trials contained one or zero predicted interlopers among 
the entire candidate sample of 45 stars - one red giant 
branch star. 90% of all trials had 7 or fewer interlopers 
(3 bRGB, 2 fRGB, and 2 MS/BHB). 

To sanity check this prediction, we also use the above 
method to predict the number of Milky Way stars ex- 
pected in our dataset outside the velocity cut for mem- 
bership in our original CM-selected sample. At 50% 
(90%) confidence, the number of stars with v < —30 
km s _1 or v > km s" 1 we predict to be in the bRGB, 
fRGB, and MS/HB regions of our spectroscopic dataset 
arc < 5, 4, and 4 (8, 7, and 7), respectively. Figure [3] 
shows that there are actually 5, 5, and 3 stars with 
V < -30 km s" 1 or v > km s" 1 in the RGB, fRGB, 
and MS/HB regions of our spectroscopic dataset. Our 
technique thus accurately predicts the number of stars 
in our spectroscopic catalog at outlying velocities. 

We also use a simpler technique that does not rely on 
the absolute calibration of the Besancon model to predict 
the number of interloper stars among the 45 candidate 
members. This "scaled histogram" approach instead re- 
lies on assuming that the 13 bRGB, fRGB, and MS/HB- 
colored stars stars at outlying velocities are all Milky 
Way foreground stars. Treating each of these three re- 
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Fig. 4. — Left panel, top: The fraction of all stars, as a function of elliptical distance from the Wil 1 center, that pass color-magnitude 
criteria for Wil 1 membership and also end up in our final spectroscopic sample with well-measured velocities. This target efficiency is 
shown for three different subpopulations of Wil 1 stars - bRGB candidates, fRGB candidates, and BHB/MS candidates. Left panel, bottom: 
The probability that any individual star is an interloper, as a function of distance. Right panel: The predicted number of Milky Way 
interlopers in the sample of 45 candidate members in 10 5 Poisson realizations of the foreground. Color coding is the same as the left panel. 
90% of all trials had 7 or fewer interlopers (3 bRGB, 2 fRGB, and 2 MS/BHB). 



gions of the CMD separately, we determine the relative 
numbers of Milky Way stars in the Besancon model with 
velocities inside and outside the velocity cut for mem- 
bership and apply this scaling factor to the numbers of 
stars we observe outside the velocity cut for membership. 
We use these fractional numbers of stars to generate 10 5 
Poisson realizations of the predicted number of bRGB, 
fRGB, and MS/HB interlopers and find that 50% (90%) 
of trials contain < 1, 1, and (3, 2, and 2) interlopers 
respectively. This technique thus yields very similar re- 
sults as the full simulations, predicting a median of only 
2 interlopers, with fewer than or equal to 7 at 90% con- 
fidence. 

We summarize these predictions for interlopers among 
the 45 candidate members in Table 4. 

3.3. Identifying Interloper Stars in Wil 1 Sample 



The calculations in S I3.2l revealed that we expect 1 — 7 
Milky Way interlopers among the 45 Wil 1 candidate 
stars identified in §3.1. When broken down by subsets 
of stellar population, we predict 1 — 3 interlopers among 
the 5 bRGB candidates, — 2 interlopers among the 10 
fRGB candidates, and — 2 interlopers among the 30 
MS/HB candidates. We now attempt to identify these 
interlopers. Because the fractional contamination of the 
candidate MS/HB Wil 1 members is small compared to 
that of RGB stars, and because we do not have a rea- 
sonable spectroscopic [Fe/H] indicator for the MS/HB 
candidates, we only look for the 1 — 5 interlopers with 
RGB magnitudes and colors. 

We use CaT reduced equivalent width, W (calculated 
using the Rutledge definition described in §2.5), to flag 
possible interlopers. All recent spectroscopic studies of 
Milky Way dwarf galaxies us e some metallicity in dicator 
to select member stars (e.g. iWalker et al.1l2009[ ). This 
selection means that the abundance spread we will infer 
for Wil 1 in §4.1 is necessarily a lower limit. We choose 
not to use velocity and position information to perform 
a likelihood analysis for member selection. We will show 



evidence that Wil 1 is both spatially and kinematically 
disturbed in § 4.2 and 5.2. We therefore do not want 
to assert that Wil 1 stars follow an exponential spatial 
distribution and Gaussian velocity distribution when dis- 
criminating member stars from interlopers. 

Figure M shows W of the 15 candidate Wil 1 RGB 
stars as a function of r$ magnitude. This figure shows 
a large spread of W, with a big gap between the more 
metal-poor and more metal-rich stars at bright magni- 
tudes. Given the gap at bright magnitudes, we hypothe- 
size that the 4 higher W' (more metal-rich) stars on this 
figure are possible foreground stars. The dotted line at 
W = 3.9 A in Figure [5] shows our adopted separation 
between possible foreground dwarfs and Wil 1 members. 
Using this W' cut, there are 4 likely foreground RGB 
stars (Stars 3, 4, 10, and 11 in Table 2). 

The fainter stars in Figure [5] do not show the same bi- 
modal distribution of W as the brighter stars. However, 
we adopt a cut at W = 3.9 A as a reasonable way to 
flag possible interlopers because i. in §3.4 we will pro- 
vide additional support for an interloper classification of 
Stars 3 and 4 and ii. we will show that this specific W' 
cut impacts neither our conclusions about the classifica- 
tion of Wil 1 nor our conclusions about the kinenratic 
properties of Wil 1. 

We test this W cut by considering in detail the two 
brightest stars flagged as interlopers. One of these 
stars, Star 4 , has a high-resolution HET spectrum from 
ISiegel et"a l. (2008). They could not obtain a consistent 
solution for Fe I and Fe II abundances under the assump- 
tion that this star was a giant. They argue it is likely to 
be a foreground dwarf, consistent with our classification. 
Star 3 has W and SDSS colors very similar to that of Star 
4 (see Tables 2 and 5). We thus consider it likely that 
Star 3 is also a fo reground star , although this statement 
is not conclusive. iMartin et al. (2007) also classified Star 
3 as an interloper. 

In addition to the four relatively high W' stars, we 
classify the star with r ~ 19 and W' ~ 3.2 (Star 6 in 
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Fig. 5. — The reduced equivalent widths (Rutledge definition) of the calcium triplet feature measured for the 15 candidate Wil 1 member 
stars brighter than the base of the red giant branch. This figure shows a large spread in W, with a gap between the more metal-rich and 
more metal-poor stars at the brightest magnitudes. Ovcrplotted is a line at 3.9 A, our subjective W' cut between probable members and 
possible interloper stars. 



Table 2) as an interloper. Figure [3] shows that Star 6 (g — 
r = 0.39) lies blueward of the main bRGB locus in g — r, 
despite its intermediate W'. For Star 6 to reasonably 
be part of Wil l's stellar population, it would need to 
be quite metal-poor (inconsistent with its intermediate 
W) to explain its blue color. Star 6 is thus a likely non- 
member, bringing the total number of stars flagged as 
contaminants on the RGB to 5 out of 15, which is at the 
upper edge of our contamination estimates. 

To facilitate the reader reaching her own conclusions 
about the classification of the five brightest Wil 1 RGB 
candidate members, the SDSS de-reddened ugriz magni- 
tudes of these stars are compiled in Table 5. To facilitate 
the reader's comparison with earlier spectros copic st udies 
of Wil 1, we also include the star I Ds used by Sic ge l et al.l 
(|2008f ) and the velocities used bv lMartin et all (j2007f )~in 
this table. 

While we cannot rule out the possibility that there 
are relatively metal-rich stars in Wil 1, we have out- 
lined a reasonable approach to flagging likely RGB non- 
members. We exclude these five stars from our primary 
analysis in the remainder of this paper. 

3.4. A close look at the membership probability of Stars 

1 and 2 

There is an abundance of circumstantial evidence for a 
member classification for both Stars 1 and 2. In their 
high-r esolution HET spectroscopic study, iSiegel et al.1 
(200$) concluded that Star 2 (unlike Star 4) is a Wil 
1 member giant. Star l's photometric and spectroscopic 
properties are very similar to that of Star 2, providing 
some evidence that they both are true member stars. The 
positions and velocities of Stars 1 and 2 are both typical 
of those of the other 43 candidate member stars: Star 1 is 
only 0.6 projected half-light radii from the center of Wil 
1 and has a heliocentric velocity of —5.4 ± 2.2 km s — 1 . 
13% of the candidate members have more positive veloc- 
ities than Star 1. Star 2 is 2.1 projected half-light radii 
from the center of Wil 1 and has a heliocentric velocity 



of — 18.5±2.2 km s _1 . 27% of the candidate members lie 
at greater distance and 22% at more negative velocities. 
(We discuss this large velocity spread and the correla- 
tion between distance and velocity in §5). We have also 
shown it to be unlikely that we have missed any fore- 
ground contaminants in the part of Wil l's RGB that 
includes Stars 1 and 2. 

We proceed with caution and now independently test 
the hypothesis that Star 1 or Star 2 may be a Milky 
Way halo star. We will show that Star l's [Fe/H] is 
— 1.73 ± 0.12, similar to th e peak of the halo's metal- 
licity distribution function (R yan fc Norrisl lT991). We 
will then argue that Wil 1 is a dwarf galaxy, or remnant 
thereof, based the [Fe/H] spread between Stars 1 and 2. 
This classification of Wil 1 thus hinges on the member- 
ship of both Stars 1 and 2. 

We use the SEGUE survey database (jYannv et al.1 
2009) to investigate whether Stars 1 and 2 are similar 
to Milky Way field stars with similar colors, magnitudes 
and velocities. We do the same for Stars 3 and 4, which 
we believe to be interlopers. We search for stars sim- 
ilar to Stars 1 and 2 by selecting SEGUE spectra of 
stars with the following restrictive set of properties: (i) 
0.51 < (g - r) < 0.55, (ii) 1.29 < (u - g) < 1.37, 
(iii) 18 < go < 19, and (iv) a radial velocity between 
—30 and 0. We search for stars similar to Stars 3 and 
4 by selecting all SEGUE spectra of stars with: (i) 
0.63 < (g - r) < 0.67, (ii) 1.67 < (u - g) < 1.74, (iii) 
18 < go < 19, and (iv) a radial velocity between —30 
and 0. We chose these photometric cuts based on the 
SDSS DR7 magnitudes of these stars (in Table 5). The 
reported SEGUE Stellar Parameter Pipeline quantities 
for the stars satisfying the Star 1/2-likc (Star 3/4-likc) 
criteria have a median [Fe/H] = —0.65 (—0.7) and log 
g = 4.42 (4.47), indicating that the samples are domi- 
nated by thick disk dwarfs, as expected for the region of 
the color-magnitude diagram we are studying. 

We downloaded the spectra of the 355 unique SEGUE 
stars satisfying Star 1/2-like criteria and of the 35 unique 
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SEGUE stars satisfying the Star 3/4-like criteria. One 
spectrum of the 355 was excluded from analysis because 
it appeared to be flawed based on a visual inspection. 
We measured the CaT EW of each SEGUE spectrum in 
the same manner as for the Wil 1 stars, using the radial 
velocities given in the SEGUE database (median error 
= 2.5 km s^ 1 ). The distribution of CaT EW for the 
SEGUE stars is reasonably well-described as a Gaussian 
with a mean of 4.1 A and a standard deviation of 0.6 A. 
The histograms are shown in Figure [5] 

No star in the Star 1/2- like SEGUE sample has a 
CaT EW < 2.2 A. The Rutledge CaT EW of Star 1 is 
1.68 ± 0.3A and Star 2 is 1.45 ± 0.3l - outliers from the 
SEGUE stars. The CaT EW of Star 1 is 2ct lower than 
the lowest of 354 stars in the SEGUE sample, showing 
that it is a true outlier from Milky Way foreground stars. 
Conversely, the stars in the Star 3/4-like SEGUE sample 
are very similar to Stars 3 and 4, which have CaT EW 
of 4.24 ± 0.3 A and 4.64 ± 0.3 A respectively. 

The spectral abundances of Stars 1 and 2 provide addi- 
tional evidence against the hypothesis that either of them 
is a field Milky Way halo star. Using the spectral synthe- 
sis method of KGS08, described in §2.5, [Fe/Hlg tor i = 
-1.73±0.12, [Fe/H] st ar2 = -2.65±0.12, [Ca/Fe] stQrl = 
-0.4±0.18, and [Ca/Fe] star2 = +0.13±0.28. Although 
its [Fe/H] is very similar to the [Fc/H] of field halo stars, 
the [Ca/Fe] of Star 1 is much lower than that of typical 
halo stars (|Venn et al.l 12004ft . The [Fe/H] of Star 2 is 
also much lower than that of typical halo stars. In § 14.11 
we will discuss the abundances of Stars 1 and 2 in more 
detail. 

Now that we have shown that neither Star 1 nor Star 
2 is likely to be a foreground star, we compare their col- 
ors and absolute m agnitudes with isochrones from the 
iDotter et al.l (|2008f ) library, and with the colors and ab- 
solute magnitudes of stars in the Draco dSph in Figures [7] 
and [5] We do this comparison to determine whether 
the very similar 5-color SDSS photometry of Stars 1 
and 2 (Table 5) can be consistent with the hypothesis 
that they are at the same distance but have very dif- 
ferent metallicities. The Draco stars (black triangles) 
overplotted with Star 1 are those with —2 .0 < [Fe/H] 
< —1 -7, as spectroscopically measured by iKirbv et al.l 
(2010). The Draco stars overplotted with Star 2 are 
those with —3.4 < [Fe/H ] < — 2.3, as spectroscopically 
measured by IKirbv et al.l (|2010( ). 

These figures show that, for the most metal-poor stars, 
the models are largely consistent with the data in g — r 
and u — g, but are too blue in g — i and u — z. This dis- 
crepancy between m odels and data i s qualitatively sim- 
ilar to that found by lAn et all (|2009f ) , who showed that 
isochrones in the SDSS photometric system are system- 
atically bluer than the colors of main sequence stars in 
relatively metal-poor globular clusters in g — i, g — z, 
u-g. 

These figures also show that stars belonging to Draco 
can have quite different spectroscopic metallicities, yet 
also have very similar colors in the SDSS filter set. Based 
on this comparison, we conclude that the similar 5-color 
SDSS photometry of Stars 1 and 2 is consistent with 
our hypothesis that they both belong to Wil 1 and have 
[Fe/H] that differ by 1 dex. 

Although our comparisons demonstrate that the sim- 



plest explanation for Stars 1 and 2 is that they belong 
to Wil 1 rather than the field halo, we cannot use this 
analysis to completely rule out the possibility that one is 
a halo star. For example, we could be unlucky and Wil 
1 may lie in a direction and at a distance that has an 
excess of stellar halo structure with abundances atypical 
relative to other lines-of-sight. However, given all of the 
evidence, we conclude that Stars 1 and 2 are both Wil 1 
members. 

3.5. A comparison with other methods of identifying 
interlopers 

Other groups have used different approaches to iden- 
tify likely Milky Way interlopers among spectroscopic 
samples of dwarf galaxy stars. The expecta t ion m axi- 
mization approach, tailored by IWalker et al.l (|2009l ) for 
the study of Milky Way dwarf galaxies, combines the 
line-of-sight velocity, Mg equivalent width, and projected 
position of each star in spectroscopic samples of Fornax, 
Carina, Sextans, and Sculptor stars to derive a probabil- 
ity that each star is a dwarf galaxy member rather than 
a contaminant star. While this approach is a powerful 
tool for studying the more luminous Milky Way galaxies, 
it is not easily applicable t o the extremely low luminosity 
satellites. First, unlike the Walker et al. (2009) datasets, 
the spectroscopic samples of the extreme satellites tend 
to include stars with a wide range of surface gravities. 
Second, we want to drop the assumptions that the veloc- 
ity distributions of stars in the extreme satellites are well 
described by Gaussians and that stars follow undisturbed 
exponential profiles out to large distances. 

Another method frequently used to identify interlopers 
is the equivalent width of the Na I line at 8190.5 A. This 
line lies in our spectral range, is sensitive to both surface 
gravity and temperature, and has previously been used to 
discriminate between red gi ant branch sta r s and red main 
sequence stars. However, iGilbert et al.1 (|2006f ) demon- 
strate that while this spectral feature is very effective for 
dwarf/ giant discrimination for V — I > 2.5, it is not ef- 
fective for stars bluer than V — I = 2.0. The red giant 
members of Wil 1 in our sample have 0.4 < V — I < 1.0, 
so we cannot use the Na I line to distinguish them from 
foreground d warfs. It is possibl e that the use of this 
diagnostic by iMartin et al.l ()2007f ) in their spectroscopic 
study of Wil 1 led to the inclusion of interlopers in their 
Wil 1 sample. For exampl e, they classify our Star 4 as a 
member giant even though lSiegel et aLl ((2008) show from 
a high-resolution spectrum that it is very likely to be a 
foreground dwarf. 

3.6. Wil 1 stars at outlying velocities? 

To investigate the possibility that we missed one or 
more member stars associated with Wil 1 (bound or un- 
bound) , we looked at the five fRGB-colored stars at out- 
lying velocities. The equivalent widths of their Na I ab- 
sorption lines are similar to those of likely Wil 1 members 
at similar apparent magnitudes. However, the surface 
gravities of faint red giant branch stars associated with 
Wil 1 are not very different than those of MW foreground 
stars in that region of the CMD. Even if this surface grav- 
ity indicator was robust for stars at these g — r colors, 
it would not provide an effective discriminant. All five 
have higher CaT equivalent widths than Wil 1 members 
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Fig. 6. — Histograms of the CaT EW measured from stellar spectra obtained from the SEGUE database. SEGUE stars were selected 
as those with color, apparent magnitude, and radial velocity extremely similar to those of Stars 1 & 2 or Stars 3 & 4 (criteria listed in 
§3.4). Stars 1 and 2 are both outliers from the 354 SEGUE stars with similar properties, showing that they are each inconsistent with 
being a foreground star. Even without the other evidence that Star 1 is a Wil 1 member star, this comparison implies the chance of Star 
1 belonging to the Milky Way foreground to be less than 1/350. Conversely, Stars 3 and 4 are both very similar to the 35 SEGUE stars 
with similar properties, supporting our hypothesis that they are foreground stars. 
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Fig. 7. — Color-magnitude diagrams showing Star 1 (red dot) compared to D raco stars (black triangles) and Dotter et al. ( 2008) isochrones. 
Draco stars with —2.0 < [Fe/H] < —1.7, as measured by Kirby et al. (20lJJ) are included in the figure. All magnitudes are de-reddened 
SDSS DR7 magnitudes, converted to absolute magnitudes, la color-magnitude measurement uncertainties are shown for each star. The 
isochrone is [Fe/H] = —1.7, [ct/Fe] = -0.2 because it was the available isochrone that was closest to the measured [Fe/H]g tar i = — 1.73±0.12, 
[Ca/Fe]g tar i = —0.4 ± 0.18. We have offset the model z-isochrones by -0.06 mag and the model i-isochrones by +0.03 mag, to put them 
in the DR7 SDSS system (A. Dotter, private communication). 



at similar apparent magnitudes. We did not attempt this 
with MS/HB-colored stars because the error bars on all 
measured parameters are too large to provide a meaning- 
ful result. Although there could possibly be one, there is 
no obvious candidate for a star associated with Wil 1 in 
the outlying tails of its velocity distribution. 

4. THE NATURE OF WILLMAN 1: A DISRUPTING DWARF 
4.1. An [Fe/H] spread 



Despite the observed star-to-star variations in some 
of the light elements, the Milky Way's globular clus- 
ters generally have not been observed to have a s ignifi- 
cant dispersion in [Fe/H] (e.g. lCarretta et a!]|2009[ ). Al- 
though spectroscopic evidence exists for modest [Fe/H] 
spreads (~ 0.1 - 0.2 dex) in a small number of 
Milky Way globular clusters (e.g. M2 2 and M54 
iDa Costa etatl 120091 : ICarretta et aJ]l2010T ). only wCen 
displ ays a large (~ 1 dex) star- to-star variation in [Fe/H] 
(e.g. iNorris fc Da Costal Il995f ). Owing to its unusual 
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Fig. 8. — The same as Figure 7, but for Star 2 (red dot) plus Draco stars (black triangles) with [Fe/H] < 
as low as -3.4, lKirbv et al.ll20TOf) . The isochrone is [Fe/H] = -2.5, [o/Fe] : 
to the measured [Fe/H] star i = -2.65 ±0.12, [Ca/Fe] star2 = +0.13 ±0.28. 



2.3 (includes stars with [Fe/H] 
+0.2 because it was the available isochrone that was closest 



abundance distribution, wCen is regarded to be the re- 
maining core of an otherwise destroyed dwarf galaxy 
(|Lee et alj H99l iBekki fc Freernan] 12001 . Unlike the 
globular clusters, all Milky Way dwarf galaxies are ob- 
served to have a significant dispersion in their stars' 
[Fe/H]. Observing such a spread is thus good evidence 
that an object self-enriched within a dark matter halo — 
that it is a galaxy. 
Assessing the evidence for a metallicity spread in Wil 

1 is challenging, because even minimal foreground con- 
tamination could lead to the erroneous conclusion that a 
large [Fe/H] spread exists. In §3.3, we applied a CaT W' 
cut to reject possible foreground stars, leaving 10 highly 
probable member stars of the initial sample of 15 possi- 
ble RGB members. This selection left only Stars 1 and 

2 of the five brightest possible member stars. We then 
performed a detailed analysis in §3.4 to ensure that Stars 
1 and 2 are true Wil 1 members. 

Using the spectral synthesis method of KGS08, de- 
scribed in §2.5, [Fe/H] star i = -1.73 ± 0.12 and 
[Fe/H] stQr2 = -2.65 ±0.12, indicating that Star 2 is 0.9 
dex more iron-poor than Star 1. To show these results 
aren't sensitive to S/N, we redo this analysis after ar- 
tificially reducing the S/N in each spectrum by a fac- 
tor of two with Gaussian random noise proportional to 
^/pixclvariance. We find [Fe/H]s tar i = — 1.77±0.12 and 
[Fe/H] S tar2 = -2.86 ± 0.20. The KGS08 technique also 
yields [Ca/Fe]stan = -0.4 ± 0.18, and [Ca/Fe] Sta r2 = 
+0.13 ± 0.28. This large spread of [Ca/Fe] abundances, 
with the more metal- rich star having the lower [Ca/Fe] , 
is consistent with a scenario where Type la supernovae 
are controlling the chemical enrichment of Wil 1 after 
its first generation of stars. With only two bright RGB 
stars, it is unlikely that these spreads sample the full 
spread of abundances in Wil 1 . 



Stars 1 and 2 also have accurate photometry in SDSS 
DR7. We use the photometric metallicity calibration of 
iLenz et al.l (|1998l ): I = -0.436u + 1.129.g - 0.119r - 
0.574i + 0.198, valid in the range 0.5 < g - r < 0.8. 
Stars with larger I are more metal-poor. Stars 1 and 2 
have I = 0.137±0.032 and 0.216±0.035 respectively. Al- 
though there is no robust conversion between this Lenz 
photometric metallicity statistic and stellar metallicity, 
we use their calculations to estimate approximate photo- 
metric metallicities from the I parameter. We fit a spline 
between I and metallicity for the full set of their calcu- 
lations shown in their Figure 8. This fit yields [M/H] 
estimates of —1.8 ± 0.3 and — 3.0^Qg for Stars 1 and 2, 
respectively. Despite the significant random and (likely) 
systematic errors in this process, these values are consis- 
tent with our spectral synthesis measurements. 

Figure [9] shows the spectra for these two stars in the 
region around the CaT along with the ratio of the spec- 
tra. The spectral syntheses used to measure [Fe/H] are 
overplotted. Visually, the individual spectra appear dif- 
ferent, and weak lines can be seen in Star l's spectrum 
that are not visible in Star 2's. The ratio of the spec- 
tra clearly shows that metal lines are stronger in Star 1. 
Given that the two stars have nearly identical luminosi- 
ties and colors, this visual comparison also demonstrates 
that Star 1 is significantly more metal- rich than Star 2. 

The scatter in W' of the eight relatively faint Wil 1 
RGB members seen in Figure [S] provides additional ev- 
idence for an internal abundance spread in this object. 
Taking the measurement uncertainties on W' into ac- 
count, these stars have a mean W' = 2.4 A and a disper- 
sion = 0.7 A. 

The method of KGS08 has been calibrated for stars 
with surface gravities lower than 3.6. Of the other eight 
probable Wil 1 RGB members, only one has log g < 3.6 
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Fig. 9. — Portions of the spectra of Stars 1 and 2 (as listed in 
Table 2) around the Ca II triplet features. The bottom spectrum 
shows the ratio of these two stars; the clear features underscore the 
true difference in the abundances of these two Wil 1 member stars. 
The spectral syntheses used to measure [Fe/H] are overplotted. 

(Star 8). We find that this star has [Fe/H] = -1.92±0.21. 
Averaging the [Fe/H] of Stars 1, 2, and 8, we find a 
mean [Fe/H] of —2.1. At face value, this is inconsistent 
with a simple linear metallicity-luminosity relation for 
M ilky Way dwarf sat ellites. For example, the relation 
of IKirbv et all (1201 ID P redicts a mean [Fe/H] = -2.7. 
Conversely. iKirbv et al.l (|2011l )'s linear relation predicts 
that the typical luminosity of a Milky Way satellite with 
[Fe/H] = -2.1 is - 6.5 x 1O 4 L . This offset between 
the observed and expected [Fe/H] of Wil 1 could result if 
the metallicity-luminosity relation flattens as the lowest 
l uminosities (as also suggested by observations of Segue 
l lSimon et al.ll20lil ). Alternatively, it could result if Wil 
1 has been stripped of a lot of its stellar mass, or if the 
mean [Fe/H] based on three stars may not be an accu- 
rate reflection of the average composition of stars in this 
system. 

For the purpose of compa rison, the [Fe/H ] deriv ed for 
Stars 1 and 2 based on the iRutledge et al.l (|1997l ) cali- 
bration of CaT W are -1.97±0.17 and -2.07±0.17, re- 
spe ctively. The [Fe/H] derive d for Stars 1 and 2 based on 
the IStarkenburg et all (2010) calibration of CaT W are 
-2. 36±0.20 and -2.64±0.29 , respectively. As expected, 
the IStarkenburg et all (120101) values a re lower metallic- 
ity, because the IRutledge et all (|1997l) [Fe/H] have been 
shown to systematically overestimate [Fe/H] lower than 
-2.0 (see KGS08, Stark enburg et al. 2010 and references 
therein). What is initially unexpected is that the [Fe/H] 
values of Stars 1 and 2 as derived using the CaT W as 
an [Fe/H] indicator underestimate the spread in [Fe/H] 
between Stars 1 and 2. This underestimate results from 
the large spread in [Ca/Fe] between Stars 1 and 2 and 
underscores another weakness of the CaT W approach 
to [Fe/H]: The CaT technique doesn't account for differ- 
ences in [Ca/Fe]. If we compare the [Ca/H] we obtain 
using the relationship provided by IStarkenburg et all 
(2010) with the [Ca/H] derived through spectral syn- 
thesis, we find excellent agreement: [Ca/H]s tar i stark 
= -2.11 ± 0.20 and [Ca,/E]stari,Kirby = -2.13 ±'0.13, 



[Ca,/B]star2, Stark — —2.39 ± 0.29 and [Ca,/B]star2.Kirby 

= -2.52 ±0.25. 

In summary, Wil 1 stars exhibit a substantial [Fe/H] 
spread. The [Fe/H] spread presented here is an under- 
estimate if we have thrown out true member stars as a 
result of our W' member selection criterion. We thus 
conclude that Wil 1 is a dwarf galaxy, or the remnants 
thereof. 

4.2. Tentative evidence for a disturbed morphology 

The spatial distribution of Wil 1 stars dis- 
plays tentative evidence for multi-direct ional features 
(Willman et all 120061: IMartin et all 120071 although see 
iMartin et al.ll2008bl and lWalsh et al.ll2008l for discussion 
of shot noi se and morphology) a nd has a moderately high 
ellipticity (jMartin et al.ll2008bl e = 0.47 ± 0.08). Addi- 
tionally, the spatial distribution of Wil l's spectroscopic 
member stars provides tenuous evidence for an extended 
spatial distribution. Figure [TU] shows that two of the 40 
Wil 1 spectroscopic members lie at ~ 5-Rhaif • Only 1% of 
the stars following an exponential distribution lie beyond 
4 Rh&it- To test the likelihood of two outlying members 
occurring by chance in a system well described by an ex- 
ponential distribution, we ran Monte Carlo simulations 
of the expe cted distribution of Wil 1 stars, assuming the 
(|Martin et al.l [2008b) structural parameters (Rh = 2.3'; 
e = 0.47) and the target efficiency in Figure 5. In only 
1.1% of simulations were there two or more stars between 
4 and 6 Rh- Using the simulations described in §3.2, we 
find it unlikely that either of these stars is a contaminant. 
For example, there is a < 2% chance that they are both 
foreground stars. These probabilities all depend on the 
assumption that Wil 1 is perfectly described by an expo- 
nential function (a larger fraction of a Plummer galaxy's 
light resides at large dista nce) and that its scal e length is 
exactly that measured by (jMartin et a l. 2008b). We thus 
consider these distant members to provide only tenuous 
evidence that Wil 1 may have an excess of stars at large 
radii, which could imply ongoing tidal stripping. 

Each individual piece of evidence for a possible dis- 
turbed morphology of Wil 1 is not remarkable. When 
combined, they provide a reasonable (but still tentative) 
basis for believing that Wil 1 has been structurally dis- 
turbed, perhaps by the tidal field of the Milky Way. 

5. THE KINEMATICS OF WIL 1 

5.1. Dynamical evidence for a high dark matter 
content? 

We begin by characterizing the global systemic veloc- 
ity and velocity dispersion of Wil 1 with a sin gle number, 
using the maximum likelihood method of IWalker et al.l 
(2006). This method assumes an intrinsic Gaussian dis- 
tribution for the line-of-sight velocities of Willman 1, and 
adds in quadrature to the intrinsic dispersion the mea- 
surement error on each star. Though it is possible that 
the system could be bound, in dynamical equilibrium 
and have a distribution function described by mild de- 
viations from Gaussianity at all radii, the fact that the 
measurement errors are independent and Gaussian makes 
the Gaussian form of the likelihood an adequate general 
description of the system. Using the sample of 40 proba- 
ble member stars, v sys = —14.1 km s _1 ± 1.0 km s _1 and 
a v = 4.0 km s _1 ± 0.8 km s . If we also include the 5 
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1 calculated in a running window of nine stars, not in- 
cluding those identified with high [Fe/H]. The running 
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Fig. 10. — The spatial distribution of the 40 Wil 1 member 
stars. For reference, dotted lines show 1, 2, 3, 4, 5 and 6-Rhalf 
around the c enter of the object, oriented at a position angle of 77° 
UMartin et al.ll2008al) . 

more metal-rich stars flagged as probable contaminants 
in §3.3, v sys = —12.8 km s _1 ± 1.0 km s _1 and a v = 4.8 
km s ± 0.8 km s . 3<r clipping removes no stars from 
either the 40 or the 45 star sample. 

Given the line-of-sight velocity dispersion as calcu- 
lated above, and assuming dynamical equili brium, it is 
straig htforward to estimate the mass of Wil 1. iWolf et all 
(2010) use the sample of 40 probable member stars pre- 
sented here and determine the mass within the half-light 
radius of of ~ 3.9t^'g x 10 5 M . This mass is relatively 
insensitive to the modeling of the velocity anisotropy pro- 
file and the parameterization of the light distribution 
of Wil 1. In terms of the central mass-to-light ratio, 
(M/L)y , of 770^440- A calculation including the 5 appar- 
ently metal-richer stars would imply a 40% larger mass. 
Even if all of the other assumptions were robust, we be- 
lieve this higher inferred mass would be erroneous for the 
reasons given in §3.3. 

5.2. An unusual and inconclusive kinematic distribution 

We now investigate the kinematic distribution of Wil 
1 stars in more detail. Is it even reasonable to char- 
acterize the systemic velocity and velocity dispersion of 
Wil 1 with a single number? Figure [TT] shows the line-of- 
sight velocities of the probable member stars (red circles) 
as a function of distance. The top panel shows elliptical 
distance from the center and the bottom panel shows dis- 
tance along the major-axis relative to the center. Blue 
squares in both panels show the 4 stars within 3.5rhaif 
with CaT W' > 3.9 A that we classified as likely in- 
terloper stars. Two member stars and one CaT W' > 
3.9 A star lie beyond the edge of both panels. We ex- 
clude these distant stars from further kinematic analysis 
because they may be remnants that have been stripped 
from the main body of the object (see §4.2). 

Stars in the top panel of Figure [TT] display an unusual 
kinematic distribution: central stars have velocities sys- 
tematically offset from those of more distant member 
stars. The solid line shows the systemic velocity of Wil 



to the outskirts of Wil 1 . The dotted line of Figure [TTJ 
shows the same as the solid line, but including the likely 
interloper stars. This line shows that our result is not 
affected by our criterion for interloper identification. It 
is tempting to invoke a tidally-disrupting scenario may 
provide an explanation of this behavior in the systemic 
velocity, however this velocity distribution might not be 
easily explained by existing models of tidally stripped 
dwarf galaxies. For example, iKlimentowski et al.l (2009) 
show that different viewing angles of a stripped dwarf 
galaxy are expected to reveal a symmetric velocity dis- 
tribution of unbound stream stars. We note that the 
details of this systemic velocity profile depends on the 
center used for Wil 1 and on the running window size. 
While we calculated Wil l's center (§2.1) as accurately 
as we could and the window size used here was not fine 
tuned to produce this result, if our center is inaccurate 
or if we choose a larger running window, we may not see 
such a striking shift in systemic velocity with distance. 
Given the relatively small kinematic sample from Wil 1 
and its extreme faintness it is difficult to conclusively 
interpret this systemic velocity variation. 

The bottom panel of Figure [TTJ shows that this un- 
usual kinematic distribution is not obviously a result of 
ordered rotation. Visually, a model with velocity that 
is a linear function of major axis position (as expected 
for ordered rotation) does not provide a good fit to these 
data. Additionally, even if the system were rotating with 
a reasonable velocity, this signal would not be able to be 
discerned with such a small sample of discrete velocities 
(lStrigarill20l7fl . 

The variation in the systemic velocity with radius 
makes it difficult to robustly define the velocity disper- 
sion of Wil 1. In fact we find the velocity dispersion 
profile of Wil 1 to be highly dependent on the running 
window or binsize used to calculate the profile. We there- 
fore choose not to show a dispersion profile here. We note 
for the 40 Wil 1 members that i) the velocity dispersion 
of the innermost 8 stars stars is equal to km s _1 with 
an uncertainty of 2.1 kin s _1 ; ii) the velocity disper- 
sion of the 9 stars with distances between 1.0 and 1.5 r c n 
is also equal to km s _1 , with an uncertainty of 2.5 
km s _1 ; and iii) the velocity dispersion of the outermost 
9 or ten stars shown in Figure 8 is equal to km s" 1 , 
with an uncertainty of 2.1 km s _1 . The details of this 
result is sensitive to the center used for Wil 1. 

We now explicitly address each of the assumptions nec- 
essary for using the kinematics to interpret a high mass- 
to-light ratio for Wil 1. While these assumptions may be 
reasonable for many of the Milky Way's companions, it 
is not yet clear that they are reasonable assumptions for 
Wil 1. 

• All 40 stars are physically associated with Wil 1 

As stated in §3, There may be a small number of 
Milky Way stars remaining in our sample of 40 
likely Wil 1 members, likely among its apparent 
MS members. The velocity dispersions of the Milky 
Way halo and thick disk are larger than that mea- 
sured for Wil 1, so stars belonging to them may 
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Fig. 11. — Heliocentric line-of-sight velocities of probable Wil 1 member stars high CaT W' stars classified as possible interlopers in 
§3.3. Top panel: Velocities as a function of 2D elliptical distance from Wil l's center, including the 38 probable members (of 40) and the 
4 possible high [Fe/H] interlopers (of 5) with r B \\ < 3.5rh a if. Overplotted is the systemic velocity of Wil 1, calculated in a 9 star running 
window, calculated both with and without the possible interloper stars. Lower panel: The line-of-sight velocities of stars as a function 
of their ID projected distance along Wil l's major axis, including the 39 probable members and the 5 possible high [Fe/H] interlopers 
with d ma j or < 3.5 80 pc. The dotted line highlights the center of Wil 1. These panels show that neither is Wil l's systemic velocity well 
described by a single value, nor is there clear evidence of rotation. 



somewhat artificially inflate the global velocity dis- 
persion measured for Wil 1. 

• All stars associated with Wil 1 are bound and in 
dynamical equilibrium 

Our sample of 40 contains too few stars to ro- 
bustly check whether the velocity distribution of 
stars in Wil 1 is consistent with dynamical equilib- 
rium. Its unusual spatial and velocity distributions 
(Figure [TT1 and discussed above), excess of spectro- 
scopic members at large distance (§3.4), and its 
relatively high [Fe/H] for its luminosity (§5) may 
indicate that Wil 1 is a disrupted or disrupting ob- 
ject. If this were the case then the assumption of 
dynamical equilibrium would not be valid. It is 
plausible that the 40 member stars reported here 
may contain some (many) unbound stars. 

• Contribution to velocities from binary stars 

Given the small velocity dispersion that we have 
measured, and the fact that the measurement un- 
certainities are similar to the measured velocity dis- 
persion, it is possible that binary stars may be in- 
flating the global line-of-sight velocity dispersion 



measured for Wil 1 stars. iMinor et al.l (|2010f ) sim- 
ulated the effect of binary stars observed veloc- 
ity dispersions of dwarf spheroidal galaxies. They 
concluded that, while binaries do inflate the ob- 
served velocity dispersion of systems such as Wil 
1, they are not expected to have inflated the ve- 
locity dispersion of a system with er f, s = 4 km/sec 
by more than 0.8 km/sec over its ^intrinsic- Al- 
ternative binary models have suggested that the 
observed dispersions of systems with intrinsic dis- 
persions of only a few tenths of a km/sec could 
be even more affected by the presence of binaries 
(A. McConnachie, private communication). Both 
modeling and repeated observations of individual 
Wil 1 member stars would be necessary to defini- 
tively conclude the effect of binaries on its observed 
velocity dispersion. 

Symmetric velocity distribution 

The distribution of line-of-sight velocities, relative 
to the mean velocity, is symmetric for any equilib- 
rium galaxy model. While noting again that our 
sample size is too small to reach any statistically- 
robust conclusion, there are initial hints of asym- 
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metry in the line-of-sight velocity distribution (See 
Fig. 8). If this asymmetric distribution persists 
with future data sets then it would support the 
hypothesis that Wil 1 is not a dynamically equili- 
brated system. 

6. DISCUSSION AND CONCLUSIONS 

The DEIMOS spectroscopic study presented here has 
revealed several new insights into the unusual Wil 1 ob- 
ject and has underscored the importance of careful fore- 
ground characterization when studying the least lumi- 
nous Milky Way companions. We have shown that the 
velocity, color, and magnitude overlap of Wil l's stellar 
population with foreground Milky Way stars make this 
object particularly difficult to study. We thus performed 
detailed Monte Carlo simulations to calculate the pos- 
sible foreground contamination. We used CaT W plus 
color information to identify a total of 5 likely interloper 
stars out of the 45 possible spectroscopic members se- 
lected in §3.1. The high interloper fraction we estimated 
for Wil l's brightest candidate RGB member stars does 
not imply that Wil l's true luminosity is less than that 
derived in past st udies. Photometri c studies of ultra- 
faint dwarfs (e.g. JMartin et all I2007t ISand et all 120091: 



010) 



IMuiioz et al.l I2010D typically use a statistical definition 
of their luminosity rather than adding up the light emit- 
ted from individual stars. 

In concert with the Besancon Galaxy model, spectra 
from the SEGUE database helped confirmed the similar- 
ity of two of the five flagged interloper stars to Milky Way 
stars with similar color, magnitude and velocity. SEGUE 
spectra also provided strong support for the presence of 
two true bright red giant branch Wil 1 members (Stars 1 
and 2 in Table 2, discussed in §3.3.1). SEGUE may thus 
provide a valuable resource for future studies that aim 
to eliminate interloper stars from spectroscopic samples. 

The mean [Fe/H] of Wil l's three confirmed RGB mem- 
bers with log g < 3.6 is —2.1, with a difference of 0.9 dex 
between the most metal-poor and metal-rich star. We 
found [Ca/Fe]stari = -0.4 ± 0.18, and [Ca/Fe] Sta r2 = 
+0.13 ± 0.28 for the brightest two stars, with the more 
Fe-rich star having the lower [Ca/Fe]. As discussed in 
§ 14.11 we interpret the large [Fe/H] difference between 
these stars to demonstrate that Wil 1 is (or once was) a 
dwarf galaxy, rather than a star cluster. With r$ ~ 18.1 
for the brighter two of these RGB members, they are 
good targets for high-resolution spectroscopic follow-up 
to investigate their detailed abundance p attern s. Be- 
cause Star 2 was observed bv lSiegel et al.l (|2008|) . Star 1 
is the top priority for follow-up. 

The kinematic distribution of Wil 1 is unlike the distri- 
bution yet seen in any of the Milky Way's satellites. Its 
inner 9 spectroscopic member stars have radial velocities 
offset by 8 km s _1 from its 29 more distant members 
(excludes the two members more distant than 3.5 fhaif)- 
Neither published models of tidally disturbed satellites 
nor ordered rotation provide an easy explanation for this 
distribution. We emphasize that the exact character of 
Wil l's systemic velocity and velocity dispersion profile 



depends sensitively on the running window and the exact 
center used. We thus use this present dataset to highlight 
the unusual nature of Wil l's kinematics rather than to 
present definitive conclusions. 

Wil l's possible disturbed morphology and tentative 
excess of spectroscopic members at large distance rela- 
tive to that of an undisturbed exponential distribution 
(§ 14.21) suggests that it may have been stripped of a sub- 
stantial fraction of its stellar component. Some models 
of the tidal evolution of dark matter dominated satel- 
lites suggest that Wil 1 should presently have a high 
M/L even if it has suff ered substantial tidal evolution 
(|Peiiarrubia et al.ll2008D . However, because its dark mass 
content cannot be well constrained given the reasons ar- 
ticulated above, appropriate caution must be taken when 
attempting to compare this object to the Milky Way's 
other satellites or attempting use it to constrain the par- 
ticle nature of dark matter. 

If Wil 1 has been severely stripped of stars and the 
line-of-sight velocities of its stars do not trace the under- 
lying gravitational potential, then why does it happen to 
lie on the Mdynamicai/iy vs Ly relationship o bserved for 
Milky Way dwarf satellites (jGeha et al.H2009t IWolf et al.l 
2010)? Perhaps this is a coincidence, or perhaps the ve- 
locities of Wil l's stars do actually trace its gravitational 
potential, despite its overall unusual kinematic distribu- 
tion. Only numerical models aimed specifically to repro- 
duce Wil l's properties may illuminate which answer is 
correct. Searches for dwarfs that can reveal the presence 
(or lack) of Willman 1-luminosity objects at a wide range 
of halo distances will be needed to know for certain the 
role of environment in shaping the luminosities of the 
tiniest Milky Way satellites. 
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TABLE 1 
Keck/DEIMOS Multi-Slitmask Observations of WillmanI 



Mask a (J2000) 
Name (h : m : s) 



<5 (J2000) PA i Gxp # of slits % useful 

(°:':") (deg) (sec) spectra 



Will.1 


10:49:23 


+51:01:20 


75 


5 x 1800 


110 


58% 


Will.2 


10:49:40 


+51:01:57 


110 


5 x 1800 


94 


45% 


Will.3 


10:49:11 


+51:02:16 


20 


3 x 1800 


92 


50% 


Will.4 


10:49:24 


+51:01:02 


-5 


3 x 1800 


127 


7% 



Note. — Right ascension, declination, position angle and total exposure time 
for each Keck/DEIMOS slitmask in Willman 1. The final two columns refer to 
the total number of slitlets on each mask and the percentage of those slitlcts for 
which a redshift was measured. Mask 4 has a low efficiency because many faint 
stars were targeted. 



TABLE 2 
Data for the 45 candidate Willman 1 members 



ID 


a (J2000) 


<5 (J2000) 


r 


(g-r) 


fhelio 


S/N 


CaT EW 




(h m s) 


(° ' ") 


(mag) 


(mag) 


( km s — 1 ) 




A 


1 


10 49 18.06 


+51 02 16 


18.13 


0.57 


-5.4 ± 2.2 


84.1 


1.63 ± 0.41 


2 


10 49 52.51 


+51 03 42 


18.14 


0.58 


-18.5 ± 2.2 


85.5 


1.41 ± 0.41 


3 a 


10 49 42.87 


+51 04 22 


18.19 


0.65 


-13.3 ± 2.2 


88.1 


4.24 ± 0.39 


l a 


10 49 12.40 


+51 05 43 


18.21 


0.63 


-7.4 ± 2.2 


72.3 


4.65 ± 0.41 


5 


10 49 07.79 


+50 56 50 


18.35 


-0.04 


-11.1 ± 2.7 


34.3 


- 


6 a 


10 49 48.53 


+51 00 32 


19.05 


0.39 


-2.1 ± 2.4 


21.6 


3.65 ± 0.46 


7 


10 49 13.13 


+51 02 32 


19.81 


-0.32 


-18.2 ± 6.8 


18.5 


- 


8 


10 49 17.41 


+51 03 25 


20.07 


0.44 


-9.4 ± 3.6 


9.4 


1.39 ± 0.46 


9 


10 49 10.11 


+51 03 00 


20.50 


0.42 


-13.1 ± 3.0 


16.3 


2.38 ± 0.40 


10 a 


10 49 15.95 


+51 02 26 


20.54 


0.45 


-8.9 ± 3.9 


9.8 


4.46 ± 0.41 


ll a 


10 49 24.97 


+51 09 23 


20.60 


0.44 


-4.9 ± 2.7 


16.8 


4.51 ± 0.39 


12 


10 49 21.14 


+51 03 29 


20.82 


0.36 


-5.6 ± 5.1 


5.6 


0.75 ± 0.44 


13 


10 48 58.10 


+51 02 53 


20.86 


0.46 


-3.8 ± 2.7 


8.3 


2.79 ± 0.40 


11 


10 49 40.82 


+51 03 39 


20.94 


0.37 


-18.9 ± 3.7 


15.3 


2.60 ± 0.40 


15 


10 49 30.94 


+51 03 40 


20.96 


0.34 


-16.4 ± 3.9 


13.2 


2.21 ± 0.41 


16 


10 49 16.75 


+51 04 03 


21.00 


0.34 


-7.6 ± 3.2 


10.6 


3.69 ± 0.44 


17 


10 50 02.85 


+51 02 32 


21.27 


0.43 


-15.8 ± 5.1 


5.4 


3.12 ± 0.41 


18 


10 50 19.29 


+51 00 12 


21.39 


0.34 


-24.2 ± 2.7 


7.5 


- 


19 


10 48 57.12 


+51 02 31 


21.42 


0.26 


-15.3 ± 3.9 


7.3 


- 


20 


10 49 28.06 


+51 01 51 


21.45 


0.22 


-21.2 ± 3.1 


3.5 


- 


21 


10 49 28.06 


+51 01 51 


21.45 


0.22 


-12.9 ± 4.6 


3.1 


- 


22 


10 49 12.55 


+51 03 05 


21.49 


0.25 


-15.9 ± 5.3 


7.0 


- 


23 


10 49 26.31 


+51 03 16 


21.49 


0.22 


-10.5 ± 6.7 


8.5 


- 


21 


10 49 25.06 


+51 02 52 


21.54 


0.25 


-13.7 ± 5.1 


1.6 


- 


25 


10 49 07.67 


+51 01 46 


21.56 


0.23 


-19.1 ± 7.2 


6.9 


- 


26 


10 48 44.41 


+51 02 06 


21.57 


0.28 


-18.4 ± 4.7 


7.3 


- 


27 


10 49 51.52 


+51 01 52 


21.58 


0.36 


-16.0 ± 3.5 


8.3 


- 


28 


10 49 34.77 


+51 04 28 


21.60 


0.27 


-18.9 ± 8.8 


7.2 


- 


29 


10 49 09.93 


+51 00 56 


21.63 


0.26 


-4.8 ± 5.0 


6.1 


- 


30 


10 49 29.73 


+51 02 58 


21.64 


0.21 


-23.5 ± 3.7 


6.1 


- 


31 


10 49 54.69 


+51 03 27 


21.66 


0.27 


-8.0 ± 4.8 


4.0 


- 


32 


10 49 14.23 


+51 01 10 


21.66 


0.23 


-10.5 ± 4.8 


7.6 


- 


33 


10 49 27.64 


+51 03 32 


21.75 


0.22 


-3.6 ± 4.8 


3.8 


- 


31 


10 49 05.17 


+51 01 42 


21.80 


0.22 


-14.8 ± 4.9 


6.6 


- 


35 


10 49 22.53 


+51 05 47 


22.04 


0.27 


-23.7 ± 9.2 


2.1 


- 


36 


10 49 11.98 


+51 03 04 


22.13 


0.28 


-19.3 ± 7.8 


3.1 


- 


37 


10 49 02.81 


+50 59 40 


22.20 


0.27 


-19.7 ± 6.5 


1.8 


- 


38 


10 49 40.69 


+51 04 19 


22.26 


0.24 


-17.1 ± 5.0 


1.1 


- 


39 


10 49 14.93 


+51 02 53 


22.44 


0.28 


-10.0 ± 4.6 


3.6 


- 


10 


10 49 13.44 


+51 02 43 


22.48 


0.36 


-5.6 ± 9.0 


3.1 


- 


41 


10 49 10.91 


+51 02 14 


22.61 


0.29 


-19.7 ± 7.4 


2.3 


- 


42 


10 49 47.62 


+51 02 03 


22.62 


0.28 


-18.1 ± 8.6 


3.5 


- 


43 


10 48 54.03 


+51 01 38 


22.63 


0.31 


-16.1 ± 5.6 


3.1 


- 


11 


10 49 40.22 


+51 03 04 


22.89 


0.29 


-9.5 ± 5.5 


2.9 


- 


45 


10 49 10.96 


+51 01 32 


22.93 


0.32 


-16.9 ± 6.4 


2.1 


- 



Note. — S/N is the median per pixel signal-to-noise for each star. Velocity error bars 
were determined from measurement overlaps as discussed in § 12.41 We supply the CaT EW 
only for stars possibly in the fRGB or bRGB populations of Wil 1. 

a Star flagged as a likely non-member by the CaT EW < 2.3 A or color criterion described 
in §3.3. 
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TABLE 3 
Data for the 52 Milky Way foreground stars in our sample 



ID a (J2000) 
(h m s) 



8 (J2000) r (g - r) v he i io 

(° ' ") (mag) (mag) (kms" 1 ) 



S/N 



46 


10 50 15.74 


+51 02 22 


15.31 


0.65 


-24.4 ± 2.2 


227.1 


47 


10 49 09.55 


+50 54 53 


15.50 


0.67 


16.7 ± 2.3 


220.2 


48 


10 49 07.46 


+51 04 06 


16.27 


0.85 


6.4 ± 2.2 


172.9 


49 


10 48 38.69 


+51 00 26 


16.95 


0.52 


-119.7 ± 2.2 


132.6 


50 


10 49 54.92 


+51 00 40 


17.38 


0.65 


-53.8 ± 2.2 


87.8 


51 


10 49 27.20 


+50 59 26 


18.57 


0.55 


155.1 ± 2.3 


43.3 


52 


10 48 45.47 


+50 58 40 


19.17 


0.57 


-171.9 ± 2.3 


37.0 


53 


10 48 47.90 


+50 57 17 


19.33 


0.40 


14.8 ± 2.3 


34.4 


54 


10 49 13.10 


+51 06 26 


19.50 


0.86 


-24.5 ± 2.6 


19.0 


55 


10 50 02.13 


+51 02 04 


19.62 


1.26 


-29.5 ± 2.2 


47.5 


56 


10 49 21.51 


+51 08 26 


20.03 


1.42 


13.2 ± 2.2 


78.3 


57 


10 50 19.02 


+50 59 18 


20.03 


0.41 


-249.3 ± 2.4 


20.3 


58 


10 50 13.69 


+51 00 06 


20.04 


1.12 


-205.0 ± 2.3 


23.3 


59 


10 49 12.90 


+51 06 17 


20.17 


0.47 


14.3 ± 2.4 


18.4 


60 


10 49 28.48 


+51 01 33 


20.27 


1.41 


-23.3 ± 2.3 


24.6 


61 


10 48 55.77 


+51 01 19 


20.32 


1.14 


-157.8 ± 2.5 


26.4 


62 


10 50 05.69 


+51 01 54 


20.33 


1.40 


-10.7 ± 2.2 


28.7 


63 


10 49 56.07 


+51 00 00 


20.61 


1.40 


6.2 ± 2.2 


24.0 


64 


10 49 22.03 


+50 59 05 


20.73 


0.32 


-36.3 ± 3.6 


8.2 


65 


10 49 20.02 


+51 04 58 


20.83 


1.35 


60.7 ± 2.3 


31.3 


66 


10 49 34.57 


+51 02 39 


20.85 


0.44 


-52.8 ± 4.8 


7.9 


67 


10 49 13.34 


+51 04 28 


20.85 


1.39 


-21.9 ± 2.3 


14.7 


68 


10 49 37.15 


+51 02 37 


20.92 


1.03 


-69.5 ± 2.4 


15.0 


69 


10 50 09.54 


+50 59 22 


21.06 


1.33 


-205.3 ± 2.4 


14.6 


70 


10 50 11.50 


+51 01 11 


21.21 


0.50 


39.9 ± 3.8 


6.5 


71 


10 49 39.43 


+51 02 03 


21.28 


1.19 


-69.1 ± 2.3 


31.3 


72 


10 49 31.01 


+51 01 33 


21.34 


0.41 


-78.9 ± 4.0 


10.3 


73 


10 49 53.75 


+51 00 43 


21.35 


1.37 


-52.2 ± 2.3 


14.9 


74 


10 49 03.74 


+51 00 38 


21.38 


1.32 


-49.0 ± 2.5 


18.1 


75 


10 49 13.88 


+51 05 12 


21.47 


1.40 


-14.7 ± 3.0 


7.9 


76 


10 49 03.41 


+51 00 48 


21.53 


0.14 


7.4 ± 5.2 


6.8 


77 


10 49 21.25 


+51 09 48 


21.79 


1.33 


-100.1 ± 2.4 


19.2 


78 


10 49 52.91 


+51 03 14 


21.80 


1.30 


-37.9 ± 3.2 


6.7 


79 


10 48 48.74 


+51 03 16 


21.82 


0.46 


65.9 ± 4.7 


6.3 


80 


10 49 07.25 


+51 02 12 


21.87 


1.44 


50.6 ± 2.4 


22.9 


81 


10 49 30.19 


+51 07 23 


21.90 


1.34 


-59.7 ± 4.3 


15.4 


82 


10 49 06.01 


+51 02 51 


22.11 


0.58 


68.1 ± 8.8 


4.2 


83 


10 49 15.36 


+51 01 05 


22.15 


1.50 


-41.4 ± 2.5 


24.1 


81 


10 49 14.03 


+51 08 44 


22.22 


1.30 


-4.9 ± 4.4 


6.4 


85 


10 49 15.46 


+51 05 51 


22.23 


1.35 


-73.4 ± 2.6 


18.3 


86 


10 49 20.72 


+51 01 41 


22.23 


0.25 


8.7 ± 7.2 


3.5 


87 


10 49 39.60 


+51 02 26 


22.25 


1.03 


87.3 ± 5.3 


16.6 


88 


10 49 24.43 


+51 09 14 


22.38 


1.40 


50.4 ± 2.6 


15.3 


89 


10 50 14.71 


+50 59 46 


22.42 


0.82 


-32.9 ± 7.1 


3.1 


90 


10 49 24.62 


+51 07 56 


22.56 


1.02 


32.9 ± 3.3 


7.9 


91 


10 48 54.87 


+51 00 07 


22.57 


1.49 


-54.4 ± 3.4 


21.2 


92 


10 49 30.07 


+51 08 19 


22.61 


1.25 


-65.2 ± 3.5 


8.0 


93 


10 50 14.66 


+51 02 16 


22.61 


1.39 


19.8 ± 6.4 


5.1 


94 


10 50 07.85 


+51 02 07 


22.77 


1.29 


-200.9 ± 5.0 


6.7 


95 


10 49 05.30 


+51 00 24 


22.89 


1.20 


124.3 ± 4.7 


12.1 


96 


10 48 57.47 


+50 57 54 


23.05 


1.15 


-10.0 ± 3.1 


9.2 


97 


10 49 38.33 


+51 01 38 


23.08 


0.44 


5.7 ± 7.8 


2.1 



Note. — S/N is the median per pixel signal-to-noisc for each star. Velocity 
error bars were determined from measurement overlaps as discussed in § 12.41 



TABLE 4 
Predicted number of MW stars in the sample of 45 candidate 

MEMBERS 



bRGB 
fRGB 

MS/BHB 



observed predicted interlopers 

simulation scaled histogram 

50% (90%) confidence 50% (90%) confidence 



5 

10 
30 



<1(3) 
<0(2) 
<0(2) 



<1(3) 

<1(2 

<0(2) 



Note. — The two methods for predicting the Milky Way contamination 
arc described in §3.2. It was calculated both with a simulation based on the 
Bcsancon Galaxy model and by scaling the histogram of the number of stars 
at velocities inconsistent with membership. 
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TABLE 5 
Additional data for bright candidate Willman 1 red giants 



ID 


Other ID a 


u h 


90 


''() 


io 


zo 


«r, Martin km S lc 


[Fe/H] d 


1 
2 
3 
1 
6 


1578 
1496 
1269 


19.99 ± 0.04 
19.91 ± 0.05 
20.53 ±0.06 
20.52 ± 0.05 
20.65 ± 0.06 


18.66 ±0.02 
18.69 ±0.02 
18.81 ±0.02 
18.84 ±0.02 
19.39 ± 0.03 


18.12 ±0.02 

18.13 ±0.02 
18.16 ±0.02 
18.19 ±0.02 
19.02 ±0.02 


17.87 ±0.02 
17.85 ± 0.02 
17.92 ± 0.02 
17.98 ± 0.02 
18.87 ±0.02 


17.72 ±0.02 
17.74 ± 0.02 
17.80 ±0.02 
17.90 ±0.02 
18.85 ±0.04 


-22.0 ± 0.6 
-13.2 ± 1.0 
-10.2 ± 1.2 


-1.73 ± 0.12 
-2.65 ± 0.12 


Note. 
a ID froi 


— Stars 1 and 6 wore neither in 
n Siegel et al. (2008). 
magnitudes arc all extinction cor 


Siegel et al. (2008) nor Martin et al. (2007). 
rectcd PSF magnitudes from SDSS Data Rcleas 








b Those 


c 7 (Abazaiian 


& Sloan Digital Skv SurvevH2009T) so 



differ slightly from the magnitudes quoted in Table 2. 
c The velocities measured by Martin et al. (2007) 
Derived using the technique described in KGS08. 



